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Abstract We explore the relationship among three coronal mass ejections (CMEs), 

observed on 28 October 2003, 7 November 2004, and 20 January 2005, the type 

II burst-associated shock waves in the corona and solar wind, as well as the ar- 

J^ rival of their related shock waves and magnetic clouds at 1 AU. Using six different 

I ^ coronal/interplanetary density models, we calculate the speeds of shocks from the 

^— ( frequency drifts observed in metric and decametric radio wave data. We compare 

C^ these speeds with the velocity of the CMEs as observed in the plane-of-the-sky 

^SJ white-light observations and calculated with a cone model for the 7 November 2004 

event. We then follow the propagation of the ejecta using Interplanetary Scintillation 

(IPS) measurements, which were available for the 7 November 2004 and 20 January 

2005 events. Finally, we calculate the travel time of the interplanetary (IP) shocks 

between the Sun and Earth and discuss the velocities obtained from the different 

O data. This study highlights the difhculties in making velocity estimates that cover 

^ the full CME propagation time. 

Ci Keywords: Coronal Mass Ejections, Initiation and Propagation, Interplanetary; 

Plasma Physics; Radio Bursts, Dynamic Spectrum, Meter- Wavelengths and 
,_^ Longer, Type II; Radio Scintillation 
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\^ 1. Introduction 

(^ In order to deduce reliable values of arrival time at Earth for CME-related ejecta, 

' it is first necessary to determine the CME departure speed from the Sun where the 
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phases of slow rise, acceleration, constant velocity, and deceleration all need to be 
considered. The problem is made more difficult by the fact that Earth-directed events 
are so-called halo or partial-halo CMEs that have a quasi-symmetric appearance 
around the solar disc when viewed in white light. Thus it is more difficult to obtain 
true earthward velocities at the departure times of these events. 

A range of tools exists for the estimation of these speeds that are mostly applicable 
from the lower corona to distances up to r -^ 10 Rq. These include the observation 
of radio bursts in the dynamic spectra (principally type II bursts) and the analysis 
of white-light images obtained with coronagraphs. For values r > 50 Rq, when the 
faster CMEs are usually decelerating, the Interplanetary Scintillation (IPS) tech- 
nique can be used. While it is possible to track the expansion of EUV or X-ray 
emitting plasmas by using appropriate imagers or by directly measuring velocities 
spectroscopically, these methods have been applied to CMEs in just a few isolated 
cases. While we expect a greater use of both imaging and spectroscopy in the future, 
through data obtained by the Hinode and STEREO missions, neither of these latter 
methods has been employed for the events discussed here. 

Since it is important to understand the strengths and limitations of these methods 
before applying them to CME speed estimation, it is necessary to briefly describe 
the underlying physics. We therefore present a short survey in the next section that 
outlines the basis for use of radio-burst dynamic spectra and white-light images for 
speed measurement. For r > 50 Rq, use of the IPS technique is most appropriate, 
and so we outline the basic features of that method. We then describe the application 
of these techniques to the events that had been selected for the Sun-Earth Connection 
Workshop (in this Topical Issue), namely those of 28 October 2003, 7 November 2004, 
and 20 January 2005, and the results obtained from the analysis. We conclude with a 
general discussion of all three halo CMEs and the uncertainties in speed estimation. 



2. Speed Measuring Techniques and Underlying Mechanisms 

Frontside "halo" CMEs are directed towards the Earth and appear as diffuse clouds 
that surround the solar disk in white-light images. Some of them can be very geoeffec- 
tive (Kim et al, 2005; Howard and Tappin, 2005). Halo CMEs are often very fast, the 
measured apparent speeds are 957 km s~ on average (Yashiro et at, 2004), which 
is about twice the mean CME speed. Speed estimation can, however, be difficult 
(Michalek, Gopalswamy, and Yashiro, 2003), as CMEs are observed in scattered 
white light on the plane of the sky, which is biased by projection effects, sideways 
expansion of the structure, and its propagation towards us. Using both white-light 
and radio observations, the "true" CME speeds have in some studies been estimated 
to be higher than the measured plane-of-the-sky speeds (Reiner et al, 2003). 

Solar type II bursts are slow-drift bursts visible in dynamic radio spectra, and 
they are generally attributed to shock-accelerated electrons (Wild and Smerd, 1972; 
Nelson and Melrose, 1985). A distinction can be made between coronal type II bursts 
observed in the decimetric-metric wavelength range and interplanetary (IP) type II 
bursts observed at decametric-hectometric (DH) wavelenghts. The driving agent of 
type II bursts has also been under discussion: shock acceleration can be created 
by blast waves {e.g., a pressure pulse without mass motions driving the wave) or 
piston-driven shocks (mass propagating at super- Alfvenic speed). There has been 
discussion on whether coronal shocks can survive into the IP space {e.g., Mann et 
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al, 2003; Knock and Cairns, 2005). IP type II bursts are usually ascribed to bow 
shocks driven ahead of a CME (Kahler, 1992), while coronal type II bursts have 
shown better temporal and spatial correlation with flare waves and ejecta {e.g., 
Klassen, Pohjolainen, and Klein, 2003; Cane and Erickson, 2005; Cliver et al, 2005, 
and references therein). 

The early work on large-scale propagating transients was based on IPS measure- 
ments of a large number of radio sources at low frequencies, applicable to relatively 
large distances from the Sun (see, e.g., Hewish, Tappin, and Capper, 1985). Present- 
day measurements at higher frequencies provide CMEs' size, speed, turbulence level, 
and mass also closer to the Sun, at distances r > 50 Rq {e.g., Manoharan et al, 
1995; Tokumaru et al, 2003). The regular monitoring of IPS on a given radio source 
over several days provides variations of the solar wind speed and density turbulence 
at a large range of heliocentric distances from the Sun. IPS observations of a grid 
of large number of radio sources on consecutive days can provide three-dimensional 
images of the heliosphere at different radii. 

2.1. CME OBSERVATIONS IN WHITE LICHT 

Since the CME structures are seen in the plane of the sky, projection effects must 
be considered. Several methods for calculating radial velocities have recently been 
proposed by, e.g., Leblanc et al. (2001), Michalek, Gopalswamy, and Yashiro (2003, 
symmetric cone model), Schwenn et al. (2005, lateral expansion speed), and Michalek 
(2006, asymmetric cone model). We selected the symmetric cone model by Michalek, 
Gopalswamy, and Yashiro (2003) for estimating the "true" speeds assuming that for 
a halo CME, propagation is with constant velocity and angular width, and the bulk 
velocity is directed radially and isotropic. If the launch location of the CME is slightly 
shifted along the cone symmetry axis by a distance r with respect to Sun centre, the 
initial and final appearance of the halo CME will be at the opposite limbs on this 
axis. Using measured plane-of-the-sky velocities at diametrically opposite limbs and 
the time difference, AT between first and last appearance, the cone model may be 
used to determine r along with i) 7, the angle between the symmetry axis and the 
plane-of-the-sky, ii) a, the opening angle of the cone, and in) V, the radial CME 
velocity. 

In view of the difficulty in estimating AT, given the relatively poor cadence of 
LASCO measurements, and since the active-region launch sites were identified for 
the events, we established r from MDI magnetogram and EIT image data, thus 
allowing an estimate of 7 from the relation cos 7 = r /Rq. Values of a and V were 
then deduced from the Equations (3) and (4) of Michalek, Gopalswamy, and Yashiro 
(2003). In fact the application of this technique was possible for only one of our 
events, that of 7 November 2004. For the event of 28 October 2003, only three 
useable LASCO C3 frames were available and the velocities at opposite limbs were 
essentially the same, making it impossible to apply the cone model. In the case of 
the event of 20 January 2005, only one useful LASCO frame was available. 

2.2. RADIO OBSERVATIONS: FROM PLASMA FREQUENCY TO SHOCK 
VELOCITY 

The general idea behind type II bursts is that they are created by a propagating 
shock: Langmuir waves are excited by electron beams produced in this shock and 
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the waves are then converted into escaping radio waves (Melrose 1980; Cairns et al. 
2003, and references therein). The emission mechanism is plasma emission near the 
fundamental and second-harmonic frequencies. 

The plasma frequency fp (Hz) at the fundamental is directly related to the electron 
density Ue (cm^ ) by 



fp = 9000V^ , (1) 

and a frequency drift towards the lower frequencies shows that the electron density 
is falling. This change is usually attributed to the burst driver moving in the solar 
atmosphere towards lower densities and larger heights. Type II bursts are therefore 
a valuable tool in determining burst (shock) velocities. 

Identifying the fundamental emission lane in dynamic radio spectra can be tricky 
(see, e.g., Robinson, 1985). Firstly not all type II bursts show emission at the funda- 
mental and second harmonic (second harmonic should be near 2/j,), and due to the 
limited observing frequencies it is possible that one or the other lane is not visible in 
the spectrum. The DH type II bursts are characterised by weaker emission than their 
metric counterparts, and the emission lanes are often split into a series of patches. 
In metric bursts the harmonic band can be stronger than the fundamental, but in 
DH bursts it is usually the opposite (see Vrsnak et al, 2001 and references therein). 
Both the fundamental and second harmonic can show band-splitting (Vrsnak et 
ai, 2001; 2002; Vrsnak, Magdalenic, and Zlobec, 2004). Sometimes rapidly drift- 
ing emission stripes called "herringbones" can be observed over a type II burst 
"backbone" (see, e.g., Cairns and Robinson, 1987; Mann and Klassen, 2005). Ground- 
based observations at decimetric-metric wavelengths are also easily affected by radio 
interference. 

After calculating the electron density from the plasma frequency at the funda- 
mental, the next step is to find a corresponding height for the emitting source. 
This is done with the help of atmospheric density models. Which model to use is a 
well-known problem, see, e.g., Robinson and Stewart (1985) and references therein. 
Figure 1 shows how density depends strongly on coronal conditions: it is important 
to know if the disturbance is propagating in a less-dense equatorial region, inside a 
dense streamer region, or in even denser coronal-loop structures. Also the turbulent 
afterfiows of a previous CME can affect the densities. The most widely used density 
models are by Newkirk (1961) and Saito (1970). In the Newkirk model, the electron 
number densities stay high at large distances from the Sun since the model is a 
hydrostatic one, see Figure 1. 

The observed density values near 1 AU are much less than in the corona (see, e.g., 
Mann et al. 1999; 2003), and the large density decreases are mainly due to solar-wind 
densities. In interplanetary (IP) space a commonly-used approximation Ue ~ R~ 
leads to fp ~ R~ , and in this case the plasma density rig (cm~ ) scales as 



'^e = ^ , (2) 

where no (cm^ ) is the plasma density near Earth at 1 AU, and Rd (in AU) is the 
distance from the Sun {e.g., Reiner et al., 2001). The value of uq ranges from -^5 
cm~ around solar minimum to ten times this value at solar maximum. 

The model of Saito (1970) departs significantly in IP space from the R^ scaling. 
A revised model by Saito, Poland, and Munro (1977) is similar to that of Saito 
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Figure 1. Electron density in different parts of the corona from eclipse photometry (image 
courtesy of S. Koutchmy, see also Koutchmy, 1994). Dark blue circles; coronal densities ac- 
cording to the Saito (1970) model, red circles: densities according to the ten-fold Saito model, 
green circles: densities according to the two-fold Newkirk (1961) model, light blue circles: 
densities according to the Leblanc, Dulk, and Bougeret (1998) model. IPM-labeled boxes 
indicate electron densities at 1.3 Kq according to the IP density model, with near-Earth 
electron densities of 5 cm"'' (solar minimum) and 50 cm~^ (solar maximum). 
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(1970) in the R = 1.2-10 Rq range, but the densities get closer to the rie '--' R^ 
relation at larger distances. Also the model by Leblanc, Dulk, and Bougeret (1998) 
gives too low densities for the active region corona (the authors recommend not to 
use it below ~1.2 Rq), but it is similar to the IP model at distances > 30 Rq. A 
"hybrid-model" by Vrsnak, Magdalenic, and Zlobec (2004) is a mixture of a five-fold 
Saito and the Leblanc, Dulk, and Bougeret models, with small modifications, and 
it can be used for connecting bursts in the corona and the IP space. Both give an 
electron density of 7 cm^ near 1 AU, which is a commonly observed value during 
activity minimum. 

Table 1. From plasma frequency to heliocentric burst driver height (-Rq) with different 
atmospheric density models 
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1 
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i) no=4.5 cm-3 at 1 AU 



We have calculated atmospheric heights for some plasma frequencies using the 
basic (one-fold) model by Saito (equatorial densities), the hybrid model by Vrsnak, 
Magdalenic, and Zlobec, a two-fold Newkirk model (approximation for active region 
corona), a ten- fold Saito model (approximation for streamer densities), the Leblanc, 
Dulk, and Bougeret model, and the IP model with an average electron density of 4.5 
cm~^ at 1 AU, see Table 1. The observed electron densities near 1 AU before the 
shock arrival for our three analysed events are listed in Table 2. 

To obtain CME velocity estimates from radio observations we make the basic 
assumption that the DH (IP) type II burst emission is formed by accelerated electrons 
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near the bow shock at the leading edge (nose) of a CME. A futher assumption is 
that the outermost bright structures in white light, projected on the plane-of-the-sky, 
represent the CME height with sufficient accuracy, although we can expect a certain 
offset between the burst driver and the bow shock (Russell and Mulligan, 2002). We 
use the observed (projected) CME heights as a constraint for the atmospheric-density 
models, and select the ones that fit best with the white-light CME observations. 
As an example for the usability of this method we refer to an event presented by 
Ciaravella et al. (2005), where the CME leading edge was identified as the shock 
front with SOHO UVCS observations on 3 March 2000, at 02:19 UT. Simuhaneous 
radio spectral observations from HiRAS show a type II burst. We calculated the 
radio source height using the fundamental emission at 40 MHz (the HiRAS dynamic 
spectrum is available at the HiRAS webpage). When the shock was observed by 
UVCS at a height of 1.7 Rq, the best-fit radio source height was 1.78 Rq, calculated 
with the "hybrid" density model. 

In determining the burst driver speed a further complication is presented by 
direction of the motion against the density gradient. If the driver is propagating 
along the radial density gradient, the observed velocity is the true shock velocity, 
but a 45°angle (0) between the directions increases the true velocity by l/cos(9, to 
about 1.4 times the observed velocity (i.e. the true distance is longer). Normally the 
direction of propagation does not differ much from the radially-decreasing density, 
so this can be taken as the upper limit for velocity correction. 

At decimetric wavelengths, a common method for determining shock speeds is 
to use the observed frequency drift rate and derive the density scale height. If we 
differentiate the expressions for plasma frequency and electron density and use the 
relation / ~ y/n, we get an expression for the shock velocity {v) 

The density scale height {H) and shock velocity (v) are expressed in km and km s~ , 
observing frequency (fundamental emission, /) in MHz, and the measured frequency 

Table 2. Estimated electron density near 1 AU, around the times of the three events 

Time no Note 

3^ 



(cm ^) 



27 Oct. 


2003 22:00 UT 


0.9 


29 Oct. 


2003 04;00 UT 


5.1 
4.0 


30 Oct. 


2003 16:00 UT 


1.0 



ACE (coronal hole prior to any ICMEs) 

ACE (right before first shock, CME launched 

on 28th) 

Geotail 

ACE (right before second shock, CME lauched 

on 29th) 



06 Nov. 2004 10:00 UT 4.5 ACE (prior to first ICME) 

09 Nov. 2004 06:00 UT 1.3 ACE (right before first shock of the second ICME) 

16 Jan. 2005 00:00 UT 4.5 ACE (small coronal hole prior to first CME) 

21 Jan. 2005 07:00 UT 3.8 ACE (right before shock) 
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drift df/dt (from a logarithmic spectral plot where it should appear as a straight 
line) in MHz s"^ 

For the local density scale height {Hl, in km) we can use a barometric isothermal 
density law for unmagnetized plasma, and the local electron density {ue) can be 
calculated from a reference density (A^e) at the base of the corona: 

/ 696 000 /-, 1\\ 
ne = iVeexp(^-^^(l--)j, (4) 

where R is the estimated heliocentric distance (observed or derived from a density 
model) in units of the solar radius Rq. Details of this method can be found in, 
e.g., Demoulin and Klein (2000). In general, the method tends to give shock speeds 
smaller than those obtained directly from the density models described above. The 
barometric density law ceases to work at large heights, where the solar-wind densities 
become a dominant factor. This density change happens at frequencies near 5 MHz, 
and for this reason we use the scale-height method only for the decimetric emission 
on 20 January 2005 in this paper. 

2.3. INTERPLANETARY SCINTILLATION OBSERVATIONS 

The propagation signatures of a CME in the space outside the LASCO field of view 
can be obtained from the remote-sensing Interplanetary Scintillation (IPS) technique 
{e.g., Manoharan et al. 2001; Tokumaru et al. 2005). The IPS method exploits the 
scattering of radiation from distant radio sources (quasars, galaxies, etc.) by the 
density irregularities in the solar wind. The normalized scintillation index (</) , where 

observed scintillation index , , 

expected average scintillation index 

can differentiate between the ambient background solar wind flow and the excessive 
level of density turbulence associated with the IP transients. The values of g close 
to unity represent the imdisturbed or background condition of the solar wind and 
values g > I and < 1, respectively, indicate the increase and decrease of density 
turbulence level in the IP medium. 

Since a propagating CME produces an excess of density turbulence by the com- 
pression of the solar wind {i.e., the sheath region) between the shock and the driver 
{i.e., the CME), the portion of the IP medium in front of the CME can be identified 
and tracked in the Sun-Earth distance with the help of scintillation images. Recent 
IPS studies have revealed the increase of the linear size of the CME with distance, 
suggesting the pressure balance is maintained between the CME cloud and the am- 
bient solar wind, in which the CME is immersed. The radial profiles of some of the 
events indicate that the internal energy of the CME supports the propagation, i.e., 
the expansion (Manoharan et al., 2000; 2001). Further, IPS data have been useful to 
quantify the force of interaction experienced by CMEs in the IP medium before their 
arrival at the near-Earth environment (Manoharan et al, 2001; Manoharan, 2006). 
In other words, the CMEs moving with speeds less than the background solar-wind 
speed are accelerated (or aided by the ambient wind). On the other hand, faster 
CMEs are decelerated due to the drag force encountered by them in the IP medium 
(Manoharan, 2006). Therefore, to assess the radial evolution of a CME before its 
arrival at 1 AU, measurements at several heliocentric distances are essential. 
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In the present study, the muhi-point IPS measurements of CMEs under investi- 
gation in the Sun-Earth space have been obtained with the Ooty Radio Telescope, 
which is operated by the Radio Astronomy Centre, Tata Institute of Fundamental 
Research, India. The description of Ooty IPS observations and the method of data 
reduction procedure have been given by Manoharan (2006) and references therein. At 
Ooty, regular monitoring of scintillation is made each day for about 700 to 900 radio 
sources. The CME events on 7-9 November 2004 and 20-21 January 2005 have been 
covered by the Ooty IPS measurements. However, due to the annual maintenance of 
the Ooty Radio Telescope, the event on 28 - 29 October 2003 has not been observed 
by the IPS method. It should be noted that the IPS measurements are sensitive to 
solar-wind structures of a CME propagating perpendicular to the line of sight to 
the radio source. An IPS image, therefore constructed using the data obtained from 
a large number of lines of sight passing through different parts of the CME, will 
provide the solar offset in the sky plane similar to the measurement obtained from 
a LASCO image. 



3. Analysis and Results from the Three Halo CME Events 

3.1. 28 OCTOBER 2003 

A fuU halo-type CME was observed at 11:30 UT in LASCO C2 (leading front at 
distance 7 Rq towards the southwest) and at 11:42 UT in C3 (leading front at 
distance 10 Rq in the same direction), see Figure 2. These two measurements give 
a CME speed of 2900 km s~^ In the LASCO CME Catalog the CME fronts are 
measured at a different angle towards the North where the heights are lower: 5.84 
Rq at 11:30 UT and 8.77 Rq at 11:42 UT. The CME speeds from the catalogue are 
2460 km s~^ (linear fit to all data points) and 2700 km s""'^ (second-order fit, speed 
near 11:30 UT). The second-order fit indicates that the CME is decelerating, and 
the speed is estimated to drop to 2000 km s~ near 30 Rq. 

A GOES X17-class flare started around 11:00 UT, superposed on an earlier event. 
Although this event is a classic halo CME, only three useable LASCO frames are 
available. For these the velocity asymmetry is not measurable and so the model of 
Michalek, Gopalswamy, and Yashiro (2003) can not be applied for the estimation of 
radial velocity. 

A bright loop front directed towards the southeast was observed before the halo 
CME, at 10:54 and 11:06 UT. The projected speed of this loop front, classified as a 
partial halo CME in the LASCO catalogue, was 1054 km s~ , which is less than half 
of the speed of the main halo event. Due to the difference in speed and direction, 
it is not clear if these two ejections are related. However, the halo event observed 
at 11:30 UT almost certainly interacted with the material from the preceding event; 
see the height-time diagram in Figure 3. 

At decimetric-metric wavelengths no clear type II emission is observed (see, e.g., 
dynamic spectra from IZMIRAN, at their webpage). At DH wavelengths. Wind 
WAVES RAD2 observed a type II burst with fundamental and second-harmonic 
emission. The burst appears in the spectrum at 11:11 UT near 14 MHz, which is 
also the upper limit for the observing frequency. Emission then drifts to the lower 
frequencies, at a rate of ~0.01 MHz s~ . At 11:16:30 UT the fundamental emission 
lane is at 4 MHz, see Figure 4, which corresponds to a plasma density of 2.0x10"'' 
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Figure 2. On the left: LASCO C2 image of the halo CME at 11:30 UT on 28 October 2003, 
with insert EIT disk image at 11:24 UT. In the middle: Corresponding running-difference 
images. The outermost CME front is located at a heliocentric distance of 7 Rq at 11:30 UT. 
On the right: LASCO C3 running-difference image at 11:42 UT, showing a symmetric halo 
structure. 
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Figure 3. GOES X-ray flux curves in the 1.5-12.5 keV and 3-25 keV energy bands, and the 
estimated heights of the halo CME on 28 October 2003. The height-time trajectory is marked 
with a solid line. An impulsive rise in soft X-rays was observed near 11:00 UT, although gradual 
rise had been recorded from 09:50 UT onwards. The flare was classified as GOES X17 class. 
Two other CMEs preceded the halo CME; their heights and times are marked with dashed 
lines. (Data from the LASCO CME Catalog.) 
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Figure 4. Wind WAVES RAD2 data at 1 - 10 MHz from 28 October 2003. Fundamental ( "F" ) 
and harmonic ("H") emission lanes of the DH type II burst are indicated in the plot. The "X" 
at 11:30 UT marks the observed distance of the CME front from the Sun centre, 7 R©, that 
corresponds to the frequency of 1.3 MHz using the Saito (1970) density model. The times and 
frequencies used for the speed calculation are indicated with dotted lines. 
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cm^ . The density model by Saito gives a corresponding atmospheric height of 3.49 
Rq, which fits quite well with the estimated CME height of 3.5 Rq at that time 
(backwards extrapolation from the observed CME leading edge heights at 11:30 and 
11:42 UT). For comparison, the hybrid, two- fold Newkirk, and ten- fold Saito models 
give unrealistically large heights, and the radio source heights given by the Leblanc 
and IP models are 0.8-2.5 Rq lower than the estimated CME front height on the 
plane of the sky. The "X" in Figure 4 at 11:30 UT marks the observed height of the 
LASCO CME front (7 Rq), corresponding to the frequency of 1.3 MHz in the Saito 
density model. 

The derived shock speeds from the radio observations using the Saito density 
model, which gives best correspondence with height, are 3630 km s~ between 
11:16:30 and 11:21:54 UT, and 2580 km s"^ between 11:21:54 and 11:30 UT, see 
Table 3. The type-II-burst lane is not fully visible at 11:30 UT, but as Figure 4 shows, 
the frequency of 1.3 MHz corresponds well to the estimated continuation of the type 
n lane. The CME speeds from the LASCO observations, approximately 2900 km s~ 
between the observations at 11:30 and 11:42 UT, and 2700 km s~ from the second- 
order fit in the LASCO CME Catalog, are not too far from the radio type H speeds 
calculated with the Saito density model. It seems evident that both the white-light 
CME and the shock driving the type II emission are decelerating. We note that a 
relatively faint streamer structure was observed towards the Southwest, which could 
affect density estimates. No interplanetary scintillation observations were available 
for this event. 

3.2. 7 NOVEMBER 2004 

The halo-CME front was first observed at a height of 5.7 Rq in the LASCO C2 
image at 16:54 UT on 7 November 2004. The front was next visible in C3 at 17:18:05 
UT at a height 9.6 Rq (LASCO CME Catalog and our analysis). The CME prop- 
agated towards the North, with a projected plane-of-the-sky speed of 1890 km s^ 
(calculated from the first two observations). The CME speeds from the catalogue are 
1760 km s~ (linear fit to all data points) and 1850 km s~ (second-order fit, speed 

Table 3. Estimated radio burst velocities on 28 October, 2003 



UT 


fp 


ne 


h Saito 


h Hybrid 


h Leblanc"* 


hlP-^ 






(cm-3) 


(Re) 


(Rq) 


(Re) 


(Rq) 


11:16:30 


4 MHz 


2.0x105 


3.49 


4.96 


2.77 


1.07 


11:21:54 


2 MHz 


4.9x10"' 


5.18 


7.07 


3.78 


2.17 


Velocity 






3630 km/s 


4530 km/s 


2195 km/s 


2365 km/s 


11:30 


1.3 MHz'' 


2.1x10"' 


7.0 








Velocity 






2580 km/s 









a) no=5.0 cm"'^ at 1 AU 

b) Frequency corresponding to the observed LASCO CME height 

Constraints: Observation of CME front at 11:30 UT: 7.0 ij0(this study), 5.84 /{©(CME 

Catalog) 

LASCO velocities: 

CME front (plane-of-the-sky) 11:30-11:42 UT: 2900 km/s 

CME front (CME Catalog, second-order fit) 11:30 UT: 2700 km/s 
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near 16:54 UT). The second-order fit indicates that the CME was decelerating at 20 
m s^^, and the speed is estimated to drop to 1600 km s^^ near 30 Rq. 

The CME was preceded by a X2.0-class flare that started at 15:42 UT, more than 
an hour before the CME was first observed. However, the first halo CME observation 
at 16:54 UT was preceded by two bad frames of LASCO C2 data, at 16:06 and 16:30 
UT. The X2.0 flare was a complex event involving at least two flares in the active 
region AR10696, and the eruption of a trans-equatorial filament. The CME appeared 
to be launched by the flare associated with the last GOES flare peak between 16:20 
and 16:54 UT (Harra et al, 2007, in this issue). 

A slow CME (225 km s^ ) was observed to propagate towards the West during 
the X-class flare, and in the C2 image at 16:54 UT in Figure 5 both CME fronts can 
be seen simultaneously. The slow CME has a height of 4.9 -Rq at that time. The 
later LASCO C3 running-difference image in Figure 5 at 17:18 UT shows the full 
halo. Figure 6 shows the height-time diagram for these structures, with the GOES 
time profile. 




Figure 5. On the left: LASCO C2 observation of tiie halo CME at 16:54 UT on 7 November 
2004, with EIT insert image at 16:48 UT. In the middle: Corresponding running-difference 
images. The outermost CME front is located northward, near 5.7 ijQ(white dashed line). The 
black dashed line on the West side of the Sun marks the location of the earlier erupted, slower 
CME front near 4.9 Rq. On the right: LASCO C3 running-difference image at 17:18 UT shows 
the asymmetric "halo" structure. 

With the availability of two LASCO C2 and six LASCO C3 frames that are 
useable, the Michalek, Gopalswamy, and Yashiro (2003) cone model (see Section 
2.1) was applied for this event. Given the identification of the launch site (active 
region that produced the X-class flare, at N09W17), the value of v/Rq = 0.287 
yields 7 = 73.3°. The values of the plane-of-the-sky velocities at the appropriate 
opposite limbs are 1146 and 660 km s~ , respectively. These, in turn, lead to values 
of a = 78.1°and V = 1488 km s^^ using Equations (3) and (4) from Michalek, 
Gopalswamy, and Yashiro. 

The 7 and a angles determine the propagation direction and angular width of 
the CME, and we can use them to deproject the true height. In the coronagraph 
images we see the projection of the CME cone's outer edge, which has an angle of 
7 — a/2 with the plane of the sky. To obtain the true heliocentric height we have 
to substract r / Rq from the projected height divided by 003(7 ~ Q!/2). For the 7 
November 2004 CME this method yields 1.2 x the projected height -0.287 Rq. If 
the projected height of the CME front at 16:54 UT was 5.7 Rq, the true height 
would have been about 6.5 Rq. Figure 5 shows that we have measured the projected 
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Figure 6. GOES X-ray flux curves in the 1.5-12.5 keV and 3-25 keV energy bands, and the 
estimated heights of the halo CME on 7 November 2004. The height-time trajectory is marked 
with a solid line. A GOES X2.0-class flare was observed to start at 15:42 UT, and the halo 
CME was first observed at 16:54 UT. A preceding CME, listed by the LASCO CME Catalog, 
is also shown (heights and times are indicated with a dashed line). LASCO C2 had bad frames 
at 16:06 and 16:30 UT, and the heights given here for the first, slower CME at 16:18 and 16:42 
UT are from LASCO C3, which blocks the lower altitudes. 



height from the slightly brightened edge towards the North. If we had made the 
measurement from the outermost bright bulk the height would be considerably less. 

The dynamic spectrum at metric wavelengths (20-70 MHz) in Figure 7 from 
GBSRBS reveals several bursts that can be connected with the bursts in the Wind 
WAVES spectrum at DH wavelengths. However, the complex event with several 
CME structures makes the identification of type H lanes difficult. Along the shock 
trajectory atmospheric densities can vary, which can explain the fast frequency drifts 
and patchy type H lanes. If the metric and DH burst emissions in this case are related, 
this could be one of the few observations of a coronal shock wave propagating into 
the IP medium. 

As the bursts in the Wind WAVES spectrum are patchy and consist of separate 
"blobs", we have selected one continuous lane of the clearly fundamental emission 
for analysis, shown in Figure 8. The LASCO CME observation at 16:54 UT is during 
this part of the type II burst, and we can compare the source heights directly 
without having to extrapolate for the CME front location. The center of the type 
II fundamental lane at 16:54 UT is at 4.5 MHz and the low frequency edge at 3.5 
MHz. The observed CME height of 5.7 Rq, at that time corresponds to a frequency 
of 1.7 MHz using the Saito density model (equatorial region), and 5.5 MHz using the 
ten- fold Saito model (high density loops or streamers). This means that the "true" 
density can be in between these two model densities, and nearer the high-density 
Saito model. An in-between model, a seven-fold Saito, would produce radio emission 
at 4.5 MHz near height 5.7 Rq. 

The derived burst-driver speed using the seven-fold Saito density model is 2955 km 
s~ . This figure is high compared to the white-light observations of the CME front 
velocity. Other density models give lower burst speeds, from 800 km s^ (Leblanc 
model) to 1765 km s~ (hybrid model), see Table 4. Using the leading edge of the 
fundamental emission lane gives slightly larger burst speeds, but the deprojected 
CME height agrees with the height given by the seven-fold Saito density model. 

The CME velocity derived from the cone model (1490 km s^^) is less than the 
plane-of-the-sky speed, but Michalek, Gopalswamy, and Yashiro do report smaller 
radial speeds in some cases. This speed agrees with the speed based on the Saito 
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Figure 7. GBSRBS dynamic spectrum at 20-70 MHz (top) and Wind WAVES spectrum at 
1.1-14 MHz (bottom) on 7 November 2004. The metric type II burst emission observed by 
GBSRBS can be interpreted as having a continuation in the DH range. 




16:30 



Figure 8. Wind WAVES RAD2 data at 1.1-12 MHz from 7 November 2004. The two lanes 
can be interpreted as fundamental and second-harmonic emission of the DH type II burst 
(indicated with dashed lines in the plot). The observed plane-of-the-sky distance of the CME 
front at 16:54 UT was 5.7 -Rq, which corresponds to a frequency of 1.7 MHz using the Saito 
(1970) density model and 5.5 MHz using the ten-fold Saito model. This frequency range is 
marked with "X"s in the plot. We made height calculations using both the center (at 4.5 
MHz at 16:54 UT) and the leading edge (at 3.5 MHz) of the emission lane at the fundamental 
(leading edge marked with a dash-dotted line). 



density model but the height does not. The best match with the height and speed 
comes from using the hybrid model. 

Figure 9 shows an IPS image of the IP medium between 18:30 UT (8 November) 
and 06:00 UT (9 November). In this "PA-heliocentric distance" image, the North is 
at the top and East is to the left. The concentric circles are of radii, 50, 100, 150, 
250 Rq. The red color code indicates the background (ambient) solar wind. The 
observing time increases from the West of the image (right side of the plot) to the 
East (left). 
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Figure 9. Ooty scintillation image (left) and CME speeds (right) from the SOHO LASCO 
CME Catalog and the measurements from the Ooty interplanetary scintillation observations 
during 7-9 November 2004. 



Table 4. Estimated radio burst velocities on 7 November, 2004 



UT 


fp 


ne 


h Hybrid 


h 7-fold Saito 


h Saito 


h Loblanc'' 






(cm-^) 


(Rq) 


(Rq) 


(Rq) 


{Re) 


16:50:19 


6 MHz 


4.4x10^ 


4.10 


4.79 


2.90 


2.33 


16:55:58 


4 MHz 


l.lxlOS 


4.96 


6.23 


3.49 


2.72 


Velocity (1 


a,ne center) 




1765 km/s 


2955 km/s 


1210 km/s 


800 km/s 


16:50:19 


4.5 MHz 


2.5x10^ 


4.70 


5.76 


3.30 


2.59 


16:54:00 


3.5 MHz 


1.5x10^ 


5.31 


6.88 


3.73 


2.87 


Velocity (1 


a,ne edge) 




1920 km/s 


3530 km/s 


1355 km/s 


900 km/s 



a) no=4.5 cm~^ at 1 AU 

Constraints: observation of CME front at 16:54 UT at 5.7 Rq, deprojected height 6.5 Rq 

LASCO and IPS velocities: 

CME front (cone model): 1490 km/s 

CME front (plane-of-the-sky) 16:54-17:18 UT: 1890 km/s 

CME front (CME Catalog, second-order fit) 18:40 UT: 1700 km/s 

CME front (IPS extrapolation) 18:40 UT: 1460 km/s 



The region of enhanced density turbulence associated with the CME is seen at 
~150 Rq. The adjacent plot shows the plane-of-the-sky speeds of the CME in the 
LASCO and IPS fields of view. It is interesting to note that the speed profile of 
the CME in the IPS field of view, when extrapolated to the LASCO field of view, 
indicates a speed of ~1460 km s~ at the heliocentric height of ~22 Rq, which 
is about the same as the speed deduced from the cone model, but about 200 km 
s~ less than the LASCO CME Catalog plane-of-the-sky speed. It is evident from 
the extended region of enhanced density turbulence seen in the IPS image that the 
CME had gone through structural changes with increasing heliocentric distance, that 
decreased the speed. 
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Figure 10. LASCO C2 observations of tiie lialo CME on 20 January 2005. At 06:30 UT (on 
tiie left) tiie CME is not yet visible. At 06:54 UT (in the middle) the CME is observed in the 
nortwest. The insert EIT images show the solar disk at 06:24 and 06:48 UT, respectively. The 
outermost CME front is located near 4.5 -Rq at 06:54 UT, measured from the centre of the 
Sun (on the right). The left and middle C2 images have reversed color scales, and the image 
on the right is a C2 running-difference image. All images suffer from intense particle hits. 



3.3. 20 JANUARY 2005 

The LASCO images for the 20 January 2005 event were hmited to only one undam- 
aged image taken at 06:54 UT, see Figure 10. Hence the cone model cannot be used 
for this event. The CME loop was directed to the Northwest, and a streamer structure 
is visible in the Southwest. The CME front was at a height of ~4.5 Rq at 06:54 UT. 
A GOES X7.1-class flare showed an impulsive rise near 06:40 UT (Figure 11), which 
caused intense particle hits and saturation effects in all instruments. Thus the speed 
of the white-light CME cannot be estimated reliably from the LASCO observations. 

In the decimetric-metric dynamic spectrum from HiRAS, several type II burst 
lanes can be observed. At 06:43 UT a patchy fundamental-second-harmonic emission 
lane pair is observed to start near 800 and 400 MHz, see spectrum at the top of 
Figure 12. At 06:54 UT a single type II lane is observed to start near 600 MHz. The 
heights of these radio sources are well below the height of the white-light CME front 
observed at 06:54 UT. Estimates for the burst heights can only be calculated with 
high-density models like the ten-fold Saito, since other density models do not have 
a solution for the high densities/high starting frequencies (Table 1). The derived 
speeds for these decimetric bursts, using the ten-fold Saito densities, are around 650 
km s~ . The density scale height method can also be used in this case. If we assume 
a base density of 10 cm~ and a heliocentric distance of 1.3 Rq for an emission 
structure at 100 MHz (in agreement with the hybrid and Newkirk density models), 
we get a local density scale height of 75 800 km. The frequency drift of the type II 
burst is around 0.6 MHz s^^ between 400 and 100 MHz (the usual metric type II 
burst drifts are between 0.1 and 1.0 MHz s~ , see Nelson and Melrose, 1985) and 
the obtained burst speed is around 940 km s^ . 

The DH type II bursts observed with the Wind WAVES instrument are faint and 
patchy. The two type-II-burst lanes, shown in the dynamic spectrum at the bottom 
of Figure 12 and indicated with dashed Hues, could both be a continuation of the 
metric emission, but due to the observational gap in the frequency range these are 
difficult to interpret. We calculated the speed for one of the burst lanes, the selected 
times and frequencies are indicated with dotted lines, and the speed varies from 750 
(Saito model) to 4690 km s^ (two-fold Newkirk model), see Table 5. These highly 
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Figure 11. GOES X-ray fiux curves in tiic 1.5-12.5 keV and 3-25 kcV energy bands, witii 
the only reliable height— time data point of the halo CME on 20 January 2005. An impulsive rise 
in soft X-rays was observed near 06:40 UT, writh a brightening in Ha. The flare wras classified 
as GOES class X7.1. 



different numbers reflect the need to have a reference height for the propagating 
shock. 

The observed LASCO CME front near 4.5 Rq at 06:54 UT cannot be used as a 
constraint here as the only LASCO observation is well before the DH type II burst 
becomes visible, and the CME height corresponds to a much lower frequency than 
where the type II appears. The frequency range corresponding to the CME height 
is marked with "X"s in the dynamic spectrum in Figure 12. This falls inside intense 
continuum emission, that would hide any type II burst. We cannot associate the 
LASCO CME front with any of the metric type II bursts either, since this would 
require extremely high densities at large atmospheric heights, which are not very 
likely. 

It is probable that if there was a DH type II burst associated with the leading 
edge of the CME, it was hidden by the intense continuum emission. Some indication 
of this is visible in the Wind WAVES RADl spectrum. A composite of the RADl 
and RAD2 observations at 80 kHz -14 MHz is shown in Figure 13, which shows 
burst patches at hectometric wavelengths between 07 and 09 UT. The times and 
heights along the indicated burst lane agree with the height of the white light CME 
at 06:54 UT (4.5 Rq) if the hybrid density model is used, as emission near 6 MHz 
comes roughly from a height of 4.1 Rq. The burst speed, calculated with the hybrid 
density model from the emission heights along the patchy lane at 3.5 MHz (07 UT) 
and 500 kHz (08 UT) is 1990 km s-\ see Table 5. 

Since the LASCO images were contaminated by the intense particle event, Gopal- 
swamy et al. (2005) combined the usable LASCO C2 image with signatures of the 
CME observed in earlier SOHO EIT data, to derive a height-time plot. Comparison 
of structures visible in the EIT difference image at 06:36 UT and in the LASCO C2 
difference image at 06:54 UT suggest a high speed, ~ 2100 km s~ . In Figure 14 
these data points are indicated by square symbols. A second-order fit to these heights 
and times would give very high speeds later on. 

Although the CME was apparently moving fast, in the IPS field of view three 
measurements could be used to track the CME at heliocentric distances >50 Rq. 
In Figure 14, a representative Ooty IPS image of the CME is shown. The estimated 
speed of the CME in the IPS field of view, on the right in Figure 14, is consistent 
with the high speed derived by Gopalswamy et al. (2005). The IPS measurements 
show that the CME decelerated from about 2500 km s~ at a solar offset of 50 Rq 
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Figure 12. HiRAS decimetric-metric spectrum (top) shows several metric type II burst lanes 
associated with the 20 January 2005 event. The Wind WAVES RAD2 dynamic spectrum 
(bottom) shows at least two possible type-II-burst lanes at decametric wavelengths. The 
type-II- burst lanes are indicated with white dashed lines in the spectral plots. The "X"s at 
06:54 UT mark the observed height of the white-light CME front at 4.5 Rq which corresponds 
to the frequency of 2.5 MHz using the Saito (1970) density model and to 7.9 MHz using the 
ten-fold Saito model. We estimated the burst speed for one of the type II lanes (the times and 
frequencies used for the calculations are indicated with dotted lines). 
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Figure 13. Wind WAVES RADl and RAD2 dynamic spectra at 06:30-09:30 UT on 20 
January 2005 (80 kHz - 14 MHz frequency range shown here) suggest that the CME-driven 
shock propagated already at hectometric wavelengths near 7 UT, see indicated burst lane 
(dash-dotted line). A clear continuous burst lane is not visible in the spectrum, but the 
individual short bursts appear at regular time intervals. The heights and times of this lane 
also agree with the white-light CME observations. 
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Figure 14. Ooty scintillation image (left) and CME speeds (right) from the estimations by 
Gopalswamy et al. (2005) and the measurements from the Ooty interplanetary scintillation 
observations during 20-21 January 2005. 



to ~1000 km s at ~1 AU. However the plot in Figure 14 suggests that the speed 
may have reached about 3000 km s^ in the interval 3 Rq to 50 Rq. 

Table 5. Estimated radio burst velocities on 20 January, 2005 



UT 


fp 


rie 


h Saito 


h 10-fold 
Saito 


h Hybrid 


h 2-fold 
Newrkirk 






(cm-3) 


(«©) 


(Re) 


(Rq) 


(Rq) 


07:16:55 


7 MHz 


6.0x10^ 


2.72 


4.86 


3.81 


5.05 


07:23:15 


5 MHz 


3.1x10^ 


3.13 


6.06 


4.46 


7.61 


Velocity 






750 km/s 


2200 km/s 


1190 km/s 


4690 km/s 


07:00 


3.5 MHz 


1.5x10^ 


3.73 




5.31 




08:00 


500 kHz 


3.1xl03 


14.62 




15.57 




Velocity 






2105 km/s 




1990 km/s 





Constraint: observation of CME front at 06:54 UT at 4.5 Rq 
LASCO and IPS velocities: 

CME front (EIT and LASCO image) estimated: >2100 km/s 
CME front (IPS) at 50 Rq: 2500 km/s 



4. Discussion on the Three Halo CME Events 

The halo CMEs on 28 October 2003 and 7 November 2004 were both preceded by a 
slower CME. Both halo CMEs were found to be decelerating after initiation. During 
the 20 January 2005 event no evident CME "interaction" was observed, and there is 
some evidence to suggest that this CME was accelerating at low atmospheric heights. 
In all three events deceleration from the initial CME speed is evident from the arrival 
times of the magnetic clouds and the measured speeds near Earth, see Table 6. The 
transit speeds of the 7 November 2004 CME/magnetic cloud are discussed in more 
detail in Culhane et al. (2007). 
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As all three CMEs were faster than the solar wind, the drag force exerted by the 
wind slowed them down considerably. Radio scintillation observations at heliocentric 
distances >50 Rq on 7-9 November 2004 and 20-21 January 2005 also record CME 
deceleration, at 3 and 30 m s~ respectively, but where this deceleration begins is 
not so clear. On 7 November 2004 the white-light CME was observed to decelerate at 
20 m s^ , which is consistent with the expectation that the CME should decelerate 
more rapidly near the Sun where the density is high, since the drag force depends on 
the density and velocity. The observed speeds near Earth from IPS measurements 
are in agreement with the measured speeds of the IP shocks and magnetic clouds. 

Wang et al. (2005) have made simulations of CME interactions, in a case where 
a slow CME is overtaken by a faster one, with the result that the faster CME is 
slowed down significantly by the blocking action of the preceding magnetic cloud. 
This implies that the travel time is dominated by the slower cloud. Therefore it is 
possible that interaction with earlier, slower CMEs affected the speeds in the first two 
events. The simulations made by Lugaz, Manchester, and Gombosi (2005) indicate 
that a CME can accelerate when it first meets the rear edge of a preceding CME 
(due to the density drop), and only decelerates later when it passes through the 
center and front of the earlier launched CME. Radio observations made at the time 
of CME acceleration would then give too high an instantaneous speed value. This 
could possibly explain the higher radio-source (shock) velocity compared to the CME 
velocity on the 7 November 2004 event, where the CME and radio source heights are 
in good agreement. 

The estimated radio source heights are also in agreement with the CME heights 
on 28 October 2003. The calculated shock velocity from the DH type II burst at 
11:16-11:21 UT (3600 km s^ ) indicates a shock speed only somewhat higher than 
the CME speed observed after 11:30 UT (2900 km s~^). The white-light CME was 
also observed to decelerate. Since the high-speed full-halo CME was propagating in 
the wake of the earlier partial-halo CME, the high momentary speed could refiect 
temporary acceleration or abrupt changes in the local density. Estimation of the 
continuation of the DH type-II-burst lane suggests that the shock was much slower 
(2600 km s^i) later on. 

The source locations at decimetric-metric wavelengths on 28 October 2003 have 
recently been analysed by Pick et al. (2005), who discovered that the radio sources 
outline an Ha-Moreton wave, which is a direct signature of a propagating shock wave 
at the photospheric level (see e.g. Uchida, 1974). The average speed of the Moreton 
wave was around 2000 km s~^. The authors also report type III bursts where the 
burst envelope moves with a speed of 2500 km s~^, in the same direction (southwest) 
as the most compact (brightest) CME structure. Therefore it is plausible that the 
burst envelope of the type III bursts and the DH type II burst were closely associated 
with the propagating CME. 

In the case of 20 January 2005 we see several type-II-burst lanes in the dynamic 
spectra at decimetric-metric and DH wavelengths, and it is possible that several 
shocks were formed during the intense event. The shock speeds derived from the 
decimetric observations (650-1000 km s~ ) are not too far from the estimated initial 
CME speed, before acceleration (as was shown in Figure 14). 
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Table 6. Estimated speeds for the magnetic 
estimates near Earth 



clouds (MCs) and the IPS velocity 



Launch 


Arrival 


Average 


MC speed 


Sun 


Notes 


time 


time 


transit 


near 


to 








speed 


Earth 


Earth 








(km/s) 


(km/s) 


(AU) 




28 Oct. 2003 


29 Oct. 


2140 


1200 


0.993 


CME launch 


11:00 UT 


05:58 UT 








to shock 
arrival 


07 Nov. 2004 


09 Nov. 


1190-1200 


800 


0.990 


CME launch 


15:53-16:20 UT 


02:00 UT 








to shock 
arrival 


07 Nov. 2004 


09 Nov. 


960-970 






CME launch 


15:53-16:20 UT 


10:00 UT 








to shock 
arrival 


07 Nov. 2004 


09 Nov. 


790-820 






to shock/ 


15:53-16:20 UT 


17:55-18:55 UT 








magnetic cloud 
arrival 


IPS near Earth 






700 






20 Jan. 2005 


21 Jan. 


1180 


675 


0.984 


flare start 


06:36 UT 


16:50 UT 








to first 
shock 


20 Jan. 2005 


21 Jan. 


1130 






flare start 


06:36 UT 


18:20 UT 








to second 
shock 


IPS near Earth 






1000 







5. Discussion on Speed Estimation 

This study of CME propagation in three different halo CME events shows that speed 
estimation can be ambiguous. We list the following items that should be taken into 
account when determining the CME propagation characteristics: 

• The plane-of-the-sky component of the speed is not simple to deproject, sup- 
posing a radial expansion from the source region, since it is biased by an 
expansion of the CME structure that is also approaching us. Furthermore, the 
CME eruption may not necessarily be radial. In addition, Thomson scattering 
is most efhcient along the perpendicular to the line-of-sight, which will lead 
to asymmetric biases in expanding structures. These biases are greatest in the 
case of frontside halo CMEs to which category all three events analysed here 
belong. Due to the biases, the C2 CME height constraints we used to deduce 
velocities using different density models may not give reliable results. The biases 
described also affect IPS measurements. The IPS distances to the Sun represent 
the orthogonal distances between the direction of the strongest disturbance and 
the Sun, which are in planes gradually turning away from the plane of the sky. 



radpap.tex; 5/02/2008; 4:32; p. 21 



Pohjolaincn ct al. 



• Simple geometric models for halo CMEs, such as the cone model, can give 
reasonable values of radial velocity provided the model constraints, e.g. constant 
velocity, are met. To solve for the two parameters - cone open angle and radial 
velocity - we need asymmetric halo CME expansion and images with good 
contrast, which reduces the number of CMEs for which the cone model may be 
used. 

• By making the assumption that shocks are formed at the nose of a propagat- 
ing CME (piston-driven bow shock with super- Alfvenic speed), we expect the 
heights of the radio shock and the heights of the white-light CME front to 
match. This might not always be the case, especially for shocks observed at 
decimetric-metric wavelengths. 

• The heights and speeds inferred from radio burst dynamic spectra need to be 
assessed critically, since they depend strongly on the selected electron density 
model. The selection of a density model includes the knowledge of certain atmo- 
spheric conditions, whether the propagation happens in equatorial quiet-Sun 
densities, in dense coronal loops or streamer regions, in undisturbed medium, 
or in the wake of earlier transients. 

• While propagation through an undisturbed interplanetary environment can be 
traced relatively well, the task becomes more challenging when propagation 
involves interaction with slower CMEs and passage through a medium that 
has been perturbed by several significant shocks in the preceding days. IP 
scintillation observations extend the range of CME observation to near Earth 
orbit, and show in particular how the faster CMEs are decelerated at distances 
>60-70 Rq from the Sun, by interaction with the solar wind plasma. 

• In-situ measurements at the LI point can determine the arrival times of shocks 
and ICME material and also allow estimates of the plasma speed. However, 
it is not always straightforward to relate particular shocks and plasma clouds 
detected near Earth to the driving CMEs that left the Sun days earlier. 
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